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Abstract. We review the use of Type la supernovae for cosmological distance deter- 
minations. Low-redshift SNe la (z < 0.1) demonstrate that the Hubble expansion is 
hnear, that Hq = 65 ± 2 (statistical) km s~^ Mpc~^, and that the properties of dust 
in other galaxies are similar to those of dust in the Milky Way. We find that the light 
curves of high-redshift (z = 0.3-1) SNe la are stretched in a manner consistent with 
the expansion of space; similarly, their spectra exhibit slower temporal evolution (by a 
factor of 1 + z) than those of nearby SNe la. The luminosity distances of our first set of 
16 high-redshift SNe la are, on average, 10-15% farther than expected in a low mass- 
density (r^A/ = 0.2) universe without a cosmological constant. Preliminary analysis 
of our second set of 9 SNe la is consistent with this. Our work supports models with 
positive cosmological constant and a current acceleration of the expansion. We address 
the main potential sources of systematic error; at present, none of them appears to 
reconcile the data with J^a — and > 0. The dynamical age of the Universe is 
estimated to be 14.2 ± 1.7 Gyr, consistent with the ages of globular star clusters. 



INTRODUCTION 

Supernovae (SNe) come in two main varieties (see reference [1] for a review). 
Those whose optical spectra exhibit hydrogen are classified as Type II, while 
hydrogen- deficient SNe are designated Type I. SNe I are further subdivided ac- 
cording to the appearance of the early-time spectrum: SNe la are characterized by 
strong absorption near 6150 A (now attributed to Si II), SNe lb lack this feature but 
instead show prominent He I lines, and SNe Ic have neither the Si II nor the He I 



lines. SNe la are believed to result from the thermonuclear disruption of carbon- 
oxygen white dwarfs, while SNe II come from core collapse in massive supergiant 
stars. The latter mechanism probably produces most SNe Ib/Ic as well, but the 
progenitor stars previously lost their outer layers of hydrogen or even helium. 

It has long been recognized that SNe la may be very useful distance indicators 
for a number of reasons; see [2,3], and references therein). (1) They are exceedingly 
himinous, with peak absolute blue magnitudes averaging —19.2 if the Hubble con- 
stant, Hq, is 65 km s~^ Mpc"^. (2) "Normal" SNe la have small dispersion among 
their peak absolute magnitudes (cr < 0.3 mag). (3) Our understanding of the pro- 
genitors and explosion mechanism of SNe la is on a reasonably firm physical basis. 

(4) Little cosmic evolution is expected in the peak luminosities of SNe la, and it 
can be modeled. This makes SNe la superior to galaxies as distance indicators. 

(5) One can perform local tests of various possible complications and evolutionary 
effects by comparing nearby SNe la in different environments. 

Research on SNe la in the 1990s has demonstrated their enormous potential as 
cosmological distance indicators. Although there are subtle effects that must indeed 
be taken into account, it appears that SNe la provide among the most accurate 
values of Hq, qq (the deceleration parameter), VIm (the matter density), and fl\ 
(the cosmological constant, Ac^/3Hq). 

There are now two major teams involved in the systematic investigation of high- 
redshift SNe la for cosmological purposes. The "Supernova Cosmology Project" 
(SCP) is led by Saul Perlmutter of the Lawrence Berkeley Laboratory, while the 
"High-Z Supernova Search Team" (HZT) is led by Brian Schmidt of the Mt. 
Stromlo and Siding Springs Observatories. One of us (A.V.F.) has worked with 
both teams, but his primary allegiance is now with the HZT. In this paper we 
present results from the HZT. 

HOMOGENEITY AND HETEROGENEITY 

The traditional way in which SNe la have been used for cosmological distance 

determinations has been to assume that they are perfect "standard candles" and 
to compare their observed peak brightness with those of SNe la in galaxies whose 
distances have been independently determined (e.g., Cepheids). The rationale is 
that SNe la exhibit relatively httle scatter in their peak blue luminosity (us ^ 0.4- 
0.5 mag; [4]), and even less if "peculiar" or highly reddened objects are eliminated 
from consideration by using a color cut. Moreover, the optical spectra of SNe la 
are usually rather homogeneous, if care is taken to compare objects at similar times 
relative to maximum brightness ( [5] and references therein) . Over 80% of all SNe la 
discovered through the early 1990s were "normal" [6]. 

From a Hubble diagram constructed with unreddened, moderately distant SNe la 
{z ^ 0.1) for which peculiar motions should be small and relative distances (as given 
by ratios of redshifts) are accurate, Vaughan et al. [7] find that 



< Ms (max) > = (-19.74 ± 0.06) + 51og {Ho/50) mag. 



(1) 



In a series of papers, Sandage et al. [8] and Saha et al. [9] combine similar relations 
with Hubble Space Telescope (HST) Cepheid distances to the host galaxies of seven 
SNe la to derive i^o = 57 ± 4 km s"^ Mpc"^. 

Over the past decade it has become clear, however, that SNe la do not constitute 
a perfectly homogeneous subclass (e.g., [1,10]). In retrospect this should have been 
obvious: the Hubble diagram for SNe la exhibits scatter larger than the photo- 
metric errors, the dispersion actually rises when reddening corrections are applied 
(under the assumption that all SNc la have uniform, very blue intrinsic colors at 
maximum; [11,12]), and there are some significant outliers whose anomalous mag- 
nitudes cannot possibly be explained by extinction alone. 

Spectroscopic and photometric peculiarities have been noted with increasing fre- 
quency in well-observed SNe la. A striking case is SN 1991T; its pre-maximum 
spectrum did not exhibit Si II or Ca II absorption lines, yet two months past max- 
imum brightness the spectrum was nearly indistinguishable from that of a classical 
SN la [13,14]. The hght curves of SN 1991T were shghtly broader than the SN la 
template curves, and the object was probably somewhat more luminous than av- 
erage at maximum. The reigning champion of well observed, peculiar SNe la is 
SN 1991bg [15-17]. At maximum brightness it was subluminous by 1.6 mag in V 
and 2.5 mag in S, its colors were intrinsically red, and its spectrum was peculiar 
(with a deep absorption trough due to Ti II). Moreover, the decline from maxi- 
mum brightness was very steep, the /-band light curve did not exhibit a secondary 
maximum like normal SNe la, and the velocity of the ejecta was unusually low. 
The photometric heterogeneity among SNe la is well demonstrated by Suntzeff [18] 
with five objects having excellent BVRI hght curves. 

COSMOLOGICAL USES: LOW REDSHIFTS 

Although SNe la can no longer be considered perfect "standard candles," they 
are still exceptionally useful for cosmological distance determinations. Excluding 
those of low luminosity (which are hard to find, especially at large distances), 
most of the nearby SNe la that had been discovered through the early 1990s were 
nearly standard ( [6], but see Li et al. [19] for recent evidence of a higher intrinsic 
peculiarity rate). Also, after many tenuous suggestions (e.g., [20-22]), convincing 
evidence has finally been found for a correlation between light-curve shape and 
luminosity. Phillips [23] achieved this by quantifying the photometric differences 
among a set of nine well-observed SNe la using a parameter, Ami5(-B), which 
measures the total drop (in B magnitudes) from maximum to i = 15 days after 
B maximum. In all cases the host galaxies of his SNe la have accurate relative 
distances from surface brightness fiuctuations or from the TuUy-Fisher relation. In 
B, the SNe la exhibit a total spread of ~ 2 mag in maximum luminosity, and the 
intrinsically bright SNe la clearly decline more slowly than dim ones. The range in 
absolute magnitude is smaller in V and /, making the correlation with Ami5(S) 
less steep than in i?, but it is present nonetheless. 



Using SNe la discovered during the Calan/Tololo survey {z < 0.1), Hamuy et 
al. [24,25] confirm and refine the PhiUips [23] correlation between Atoi5(S) and 

Mmax{B, V): it is not as steep as had been claimed. Apparently the slope is steep 
only at low luminosities; thus, objects such as SN 1991bg skew the slope of the 
best-fitting single straight line. Hamuy et al. reduce the scatter in the Hubble 
diagram of normal, unreddened SNe la to only 0.17 mag in B and 0.14 mag in V; 
see also [26]. 

In a similar effort, Riess, Press, & Kirshner [27] show that the luminosity of 
SNe la correlates with the detailed shape of the overall light curve. They form 
a "training set" of light-curve shapes from 9 well-observed SNe la having known 
relative distances, including very peculiar objects (e.g., SN 1991bg). When the 
light curves of an independent sample of 13 SNe la (the Calan/Tololo survey) are 
analyzed with this set of basis vectors, the dispersion in the V-band Hubble diagram 
drops from 0.50 to 0.21 mag, and the Hubble constant rises from 53 ± 11 to 67 ± 7 
km Mpc~^, comparable to the conclusions of Hamuy et al. [24,25]. About half 
of the rise in Hq results from a change in the position of the "ridge line" defining 
the linear Hubble relation, and half is from a correction to the luminosity of some 
of the local calibrators which appear to be unusually luminous (e.g., SN 1972E). 

By using light-curve shapes measured through several different filters, Riess, 
Press, & Kirshner [28] extend their analysis and objectively eliminate the effects of 
interstellar extinction: a SN la that has an unusually red B — V color at maximum 
brightness is assumed to be intrinsically subluminous if its light curves rise and 
decline quickly, or of normal luminosity but significantly reddened if its light curves 
rise and decline slowly. With a set of 20 SNe la consisting of the Calan/Tololo 
sample and their own objects, they show that the dispersion decreases from 0.52 
mag to 0.12 mag after application of this "multi-color light curve shape" (MLCS) 
method. The results from a recent, expanded set of nearly 50 SNe la indicate 
that the dispersion decreases from 0.44 mag to 0.15 mag (Figure 1). The resulting 
Hubble constant is 65 ± 2 (statistical) km s^^ Mpc~^, with an additional system- 
atic and zero-point uncertainty of ±5 km s^^ Mpc~^. Riess et al. [28] also show 
that the Hubble flow is remarkably linear; indeed, SNe la now constitute the best 
evidence for linearity. Finally, they argue that the dust affecting SNe la is not of 
circumstellar origin, and show quantitatively that the extinction curve in external 
galaxies typically docs not differ from that in the Milky Way (cf. [2], but sec [29]). 

The advantage of systematically correcting the luminosities of SNe la at high 
redshifts rather than trying to isolate "normal" ones seems clear in view of evidence 
that the luminosity of SNe la may be a function of stellar population. If the most 
luminous SNe la occur in young stellar populations [24,31,32], then we might expect 
the mean peak luminosity of high-rcdshift SNe la to differ from that of a local 
sample. Alternatively, the use of Cepheids (Population I objects) to calibrate local 
SNe la can lead to a zero point that is too luminous. On the other hand, as long 
as the physics of SNe la is essentially the same in young stellar populations locally 
and at high redshift, we should be able to adopt the luminosity correction methods 
(photometric and spectroscopic) found from detailed studies of low-redshift SNe la. 
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Figure 1: Hubble diagrams for SNe la [30] with velocities (km s~^) in the COBE 
rest frame on the Cephcid distance scale. Top: The objects are assumed to be 
standard candles and there is no correction for extinction; the result is o" = 0.42 
mag and i^o = 58 ± 8 km s~^ Mpc~^. Bottom: The same objects, after correction 
for reddening and intrinsic differences in luminosity. The result is a — 0.15 mag 
and ifo = 65 ± 2 (statistical) km Mpc~^. 

Large numbers of nearby SNe la are now being found by the Lick Observatory 
Supernova Search (LOSS) conducted with the 0.76-m Katzman Automatic Imaging 
Telescope (KAIT; [19,34]). CCD images are taken of about 1000 galaxies per 
night and compared with KAIT "template images" obtained earlier; the templates 
are automatically subtracted from the new images and analyzed with computer 
software. The system reobserves the best candidates the same night, to eliminate 
star-hke cosmic rays, asteroids, and other sources of false alarms. The next day, 
undergraduate students at UC Berkeley examine all candidates, including weak 
ones, and they glance at all subtracted images to locate SNe that might be near 
bright, poorly subtracted stars or galactic nuclei. LOSS discovered 20 SNe (of all 
types) in 1998 and 40 SNe in 1999, making it the world's most successful search for 
nearby SNe. The most important objects were photometrically monitored through 
BVRI (and sometimes U) filters (e.g., [33]), and unfiltered follow-up observations 
were made of all of them during the course of the SN search. This growing sample of 



well-observed SNe la should allow us to more precisely calibrate the MLCS method, 
as well as to look for correlations between the observed properties of the SNe and 
their environment (Hubble type of host galaxy, metallicity, stellar population, etc.). 

COSMOLOGICAL USES: HIGH REDSHIFTS 

These same techniques can be apphed to construct a Hubble diagram with high- 
redshift SNe la, from which the value of go can be determined. With enough objects 
spanning a range of redshifts, we can determine Qm and fl^ independently (e.g., 
[35]). Contours of peak apparent i?-band magnitude for SNe la at two redshifts 
have different slopes in the Qm^^a plane, and the regions of intersection provide 
the answers we seek. 

The Search 

Based on the pioneering work of Norgaard-Nielsen et al. [36], whose goal was 
to find SNe in moderate-redshift clusters of galaxies, the SCP [37] and our HZT 
[38] devised a strategy that almost guarantees the discovery of many faint, distant 
SNe la on demand, during a predetermined set of nights. This "batch" approach to 
studying distant SNe allows follow-up spectroscopy and photometry to be scheduled 
in advance, resulting in a systematic study not possible with random discoveries. 
Most of the searched fields are equatorial, permitting follow-up from both hemi- 
spheres. 

Our approach is simple in principle; see [38] for details, and for a description of 
our first high-redshift SN la (SN 1995K). Pairs of first-epoch images are obtained 
with the CTIO or CFHT 4-m telescopes and wide-angle imaging cameras during the 
nights just after new moon, followed by second-epoch images 3-4 weeks later. (Pairs 
of images permit removal of cosmic rays, asteroids, and distant Kuiper-belt objects.) 
These are compared immediately using well-tested software, and new SN candidates 
are identified in the second-epoch images (Figure 2). Spectra are obtained as soon 
as possible after discovery to verify that the objects are SNe la and determine their 
redshifts. Each team has already found over 80 SNe in concentrated batches, as 
reported in numerous lAU Circulars (e.g., [39], 11 SNe with 0.16 ^ z ^ 0.65; [40], 
17 SNe with 0.09 <z< 0.84). 

Intensive photometry of the SNe la commences within a few days after pro- 
curement of the second-epoch images; it is continued throughout the ensuing and 
subsequent dark runs. In a few cases HST images are obtained. As expected, most 
of the discoveries are on the rise or near maximum brightness. When possible, the 
SNe are observed in filters which closely match the redshifted B and V bands; this 
way, the X-corrections become only a second-order effect [41]. Custom-designed 
filters for redshifts centered on 0.35 and 0.45 are used by our HZT [38], when ap- 
propriate. We try to obtain excellent multi-color light curves, so that reddening 
and luminosity corrections can be applied [28,31,25]. 




Figure 2: Discovery image of SN 1997cj (28 April 1997), along with the template 
image (7 April 1997) and an HST image obtained subsequently. The net (sub- 
tracted) image is also shown. 

Although SNe in the magnitude range 22-22.5 can sometimes be spectroscopi- 
cally confirmed with 4-m class telescopes, the signal-to-noise ratios are low, even 
after several hours of integration. Certainly Keck is required for the fainter objects 
(mag 22.5-24.5). With Keck, not only can we rapidly confirm a large number of 
candidate SNc, but wc can search for peculiarities in the spectra that might in- 
dicate evolution of SNe la with redshift. Moreover, high-quality spectra allow us 
to measure the age of a supernova: we have developed a method for automati- 
cally comparing the spectrum of a SN la with a library of spectra corresponding 
to many different epochs in the development of SNe la [5] . Our technique also has 
great practical utility at the telescope: we can determine the age of a SN "on the 
fly," within half an hour after obtaining its spectrum. This allows us to rapidly de- 
cide which SNe are best for subsequent photometric follow-up, and we immediately 
alert our collaborators on other telescopes. 



Results 

First, we note that the light curves of high-redshift SNe la are broader than those 
of nearby SNe la; the initial indications [42,43] are amply confirmed with our larger 
samples. Quantitatively, the amount by which the light curves are "stretched" is 
consistent with a factor oi 1 + z, as expected if redshifts are produced by the 
expansion of space rather than by "tired light." We were also able to demonstrate 
this spectroscopically at the 2a confidence level for a single object: the spectrum 



of SN 1996bj {z — 0.57) evolved more slowly than those of nearby SNe la, by a 
factor consistent with 1 + z [5] . More recently, we have used observations of SN 
1997cx {z = 0.36) at three epochs to conclusively verify the effects of time dilation: 
temporal changes in the spectra are slower than those of nearby SNe la by roughly 
the expected factor of 1.36 [44]. 

Following our Spring 1997 campaign, in which we found a SN with z — 0.97 
(SN 1997ck), and for which we obtained //ST follow-up images of three SNe, we 
published our first substantial results concerning the density of the Universe [45]: 
= 0.35 ± 0.3 under the assumption that fltotai = 1, or Qm = —0.1 ± 0.5 
under the assumption that Qa = 0. Our independent analysis of 10 SNe la using 
the "snapshot" distance method (with which conclusions are drawn from sparsely 
observed SNe la) gives quantitatively similar conclusions [46]. 
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Figure 3 (left): The upper panel shows the Hubble diagram for the low-redshift and 
high-redshift SN la samples with distances measured from the MLCS method; see 
[48]. Overplotted are three world models: "low" and "high" Qm with Qa = 0, and 
the best fit for a fiat universe {Qm = 0.28, Q/<^ = 0.72). The bottom panel shows the 
difference between data and models from the Qm = 0.20, Qa = prediction. Except 
for SN 1997ck {open symbol; z = 0.97), which lacks spectroscopic confirmation and 
was excluded from the fit, only the 9 best-observed high-redshift SNe la are shown. 
The average difference between the data and the Qm — 0.20, JIa = prediction is 
0.25 mag. 

Figure 4 (right): Joint confidence intervals for (I^Mj^a) from SNe la [48]. The 
solid contours are results from the MLCS method applied to 10 well-observed SN la 
fight curves, together with the snapshot method [46] appfied to 6 incomplete SN la 
light curves. The dotted contours are for the same objects excluding SN 1997ck 
{z = 0.97). Regions representing specific cosmological scenarios are illustrated. 



Our next results, obtained from a total of 16 high- 2; SNe la, were announced at 
a conference in February 1998 [47] and formally pubhshed by Riess et al. [48] in 

September 1998. The Hubble diagram (from a refined version of the MLCS method 
[48]) for the 10 best-observed high- 2; SNe la is given in Figure 3, while Figure 4 
illustrates the derived confidence contours in the Qm^^a plane. We confirm our 
previous suggestion that is low. Even more exciting, however, is our conclusion 
that f^A is nonzero at the 3a statistical confidence level. With the MLCS method 
applied to the full set of 16 SNe la, our formal results are Qm = 0.24 ± 0.10 
if fitotai = 1, or flM — —0.35 ± 0.18 (unphysical) if = 0. If we demand that 
VIm = 0.2, then the best value for Ha is 0.66±0.21. These conclusions do not change 
significantly if only the 9 best-observed SNe la are used (Figure 3; Um = 0.28 ±0.10 
if fitotai = 1)- The Ami5{B) method yields similar results; if anything, the case for 
a positive cosmological constant strengthens. (For brevity, in this paper we won't 
quote the Ami5(i?) numbers; see [48] for details.) From an essentially independent 
set of 42 high- 2; SNe la (only 2 objects in common), the SCP obtains almost identical 
results [49]. This suggests that neither team has made a large, simple blunder! 

Recently, we have calibrated an additional sample of 9 high-2; SNe la, includ- 
ing several observed with HST. Preliminary analysis suggests that the new data 
are entirely consistent with the old results, thereby strengthening their statisti- 
cal significance. Figure 5 shows the tentative Hubble diagram; full details will be 
published elsewhere. 

Though not drawn in Figure 4, the expected confidence contours from measure- 
ments of the angular scale of the first acoustic peak of the cosmic microwave back- 
ground radiation (CMBR) are nearly perpendicular to those provided by SNe la 
(e.g., [50,51]); thus, the two techniques provide complementary information. A 
stunning result was already available by mid- 1998 from existing measurements 
[52,55]: our analysis of the data in [48] demonstrates that Qm + = 0.94 ± 0.26, 
when the SN and CMBR constraints are combined [54] (see also [55,56], and oth- 
ers). As shown in Figure 6, the confidence contours are nearly circular, instead 
of highly eccentric ellipses as in Figure 4. Just a few days before the Second 
Tropical Workshop, the more accurate and precise results of the BOOMERANG 
collaboration were announced [57], and shortly thereafter the MAXIMA collabora- 
tion distributed their very similar findings [58,59]; the TOCO measurements [60] 
are also relevant. The bottom line is that we appear to live in a flat universe: 
^totai = 1- Combined with the supernova results, the evidence for nonzero Q\ is 
strong. Making the argument even more compelling is the fact that various studies 
of clusters of galaxies (see summary by [61]) show that Qm ~ 0.3, so if the CMBR 
results are correct, one is led to the independent conclusion that Q\ > 0. We 
eagerly look forward to future CMBR measurements of even greater precision. 

The dynamical age of the Universe can be calculated from the cosmological pa- 
rameters. In an empty Universe with no cosmological constant, the dynamical age 
is simply the "Hubble time" (i.e., the inverse of the Hubble constant); there is no 
deceleration. SNe la yield Hq = 65±2 km s~^ Mpc~^ (statistical uncertainty only), 
and a Hubble time of 15.1 ± 0.5 Gyr. For a more complex cosmology, integrating 



the velocity of the expansion from the current epoch {z — 0) to the beginning 
{z — oo) yields an expression for the dynamical age. As shown in detail by Riess 

et al. [48], wc obtain a value of 1A.2^q^ Gyr using the likely range for {Qm,^a) 
that wc measure. (The precision is so high because our experiment is sensitive 
to roughly the difference between VIm and fl\, and the dynamical age also varies 
in approximately this way.) Including the systematic uncertainty of the Cepheid 
distance scale, which may be up to 10%, a reasonable estimate of the dynamical 
age is 14.2 ± 1.7 Gyr. 




Figure 5 (left): As in Figure 3, the upper panel shows the Hubble diagram for 
the low-^ and high- 2; SN la samples. Here, we include preliminary analysis of 9 
additional SNe la (open squares). The bottom panel shows the difference between 
data and models from the Q,m — 0.20, = prediction. 

Figure 6 (right): The HZT's combined constraints from SNe la (Figure 3) and 
the position of the first acoustic peak of the CMB angular power spectrum, based 
on data available in mid-1998; see [54]. The contours mark the 68%, 95.4%, and 
99.7% enclosed probability regions. Sohd curves correspond to results from the 
MLGS method, while dotted ones are from the Ami5(S) method; all 16 SNe la in 
[48] were used. 

This result is consistent with ages determined from various other techniques such 
as the cooling of white dwarfs (Galactic disk > 9.5 Gyr [62]), radioactive dating of 
stars via the thorium and europium abundances (15.2 ± 3.7 Gyr [63]), and studies 
of globular clusters (10-15 Gyr, depending on whether Hipparcos parallaxes of 
Cepheids are adopted [64,65]). Evidently, there is no longer a problem that the age 
of the oldest stars seems greater than the dynamical age of the Universe. 



DISCUSSION 



High-redshift SNe la are observed to he dimmer than expected in an empty Uni- 
verse (i.e., Qm = 0) with no cosmological constant. A cosmological explanation 
for this observation is that a positive vacuum energy density accelerates the expan- 
sion. Mass density in the Universe exacerbates this problem, requiring even more 
vacuum energy. For a Universe with Qm = 0.2, the average MLCS distance moduli 
of the well-observed SNe are 0.25 mag larger (i.e., 12.5% greater distances) than 
the prediction from Q\ = 0. The average MLCS distance moduli are still 0.18 mag 
bigger than required for a 68.3% {la) consistency in a universe with = 0.2 and 
without a cosmological constant. The derived value of Qq is —0.75 ± 0.32, implying 
that the expansion of the Universe is accelerating. If Qj>^ really is constant, then at 
least the region of the Universe wc have observed {z ^ 0.8) will expand eternally. 
Under the simplifying assumption of global homogeneity and isotropy, the entire 
Universe will behave in this manner. 

A very important point is that the dispersion in the peak luminosities of SNe la 
((T = 0.15 mag) is low after application of the MLCS method of [28,48]. With 
16 SNe la, our effective uncertainty is 0.15/4 ?a 0.04 mag, less than the expected 
difference of 0.25 mag between universes with = and 0.76 (and low Om); see 
Figure 3. Systematic uncertainties of even 0.05 mag (e.g., in the extinction) are 
significant, and at 0.1 mag they dominate any decrease in statistical uncertainty 
gained with a larger sample of SNe la. Thus, our conclusions with only 16 SNe la 
are already limited by systematic uncertainties, not by statistical uncertainties — 
but of course the 9 new objects further strengthen our case. 

Here we explore the major possible systematic effects that could invalidate our 
results. Of those that can be quantified at the present time, none appears to 
reconcile the data with Q\ = 0, though further work is necessary to verify this. 
Additional details can be found in [38] and especially [48]. 

Evolution 

Perhaps the most obvious possible culprit is evolution of SNe la over cosmic time, 

due to changes in metallicity, progenitor mass, or some other factor. If the peak 
luminosity of SNe la were lower at high rcdshift, then the case for Q\ > would 
weaken. Conversely, if the distant explosions are more powerful, then the case for 
acceleration strengthens. Theorists are not yet sure what the sign of the effect will 
be, if it's present at all; different assumptions lead to different conclusions [66-70]. 

Of course, it is very difficult to obtain an independent measure of the peak 
luminosity of high-;^ SNe la, and hence to directly test for luminosity evolution. 
However, we can more easily determine whether other observable properties of low- 
z and high- 2; SNe la differ. If they are all the same, it is more probable that the 
peak luminosity is constant as well — but if they differ, then the peak luminosity 
might also be affected (e.g., [66]). Drell et al. [71], for example, argue that there are 



reasons to suspect evolution, because the average properties of existing samples of 
high-z and low- 2; SNe la seem to differ (e.g., the high- 2; SNe la are more uniform). 

The local sample of SNe la displays a weak correlation between light-curve shape 
(or luminosity) and host galaxy type, in the sense that the most luminous SNe la 
with the broadest light curves only occur in late-type galaxies. Both early-type 
and late-type galaxies provide hosts for dimmer SNe la with narrower light curves 
[31]. The mean luminosity difference for SNe la in late-type and early-type galaxies 
is ~ 0.3 mag. In addition, the SN la rate per unit luminosity is almost twice as 
high in late- type galaxies as in early-type galaxies at the present epoch [72]. These 
results may indicate an evolution of SNe la with progenitor age. Possibly relevant 
physical parameters are the mass, metallicity, and C/0 ratio of the progenitor [66]. 

We expect that the relation between light-curve shape and luminosity that applies 
to the range of stellar populations and progenitor ages encountered in the late-type 
and early-type hosts in our nearby sample should also be applicable to the range 
we encounter in our distant sample. In fact, the range of age for SN la progenitors 
in the nearby sample is likely to be larger than the change in mean progenitor age 
over the 4-6 Gyr lookback time to the high-^ sample. Thus, to first order at least, 
our local sample should correct our distances for progenitor or age effects. 

We can place empirical constraints on the effect that a change in the progenitor 
age would have on our SN la distances by comparing subsamples of low-redshift 
SNe la believed to arise from old and young progenitors. In the nearby sample, the 
mean difference between the distances for the early-type (8 SNe la) and late-type 
hosts (19 SNe la), at a given redshift, is 0.04 ± 0.07 mag from the MLCS method. 
This difference is consistent with zero. Even if the SN la progenitors evolved from 
one population at low redshift to the other at high redshift, we still would not 
explain the surplus in mean distance of 0.25 mag over the Q\ = prediction. 

Moreover, it is reassuring that initial comparisons of high-redshift SN la spec- 
tra appear remarkably similar to those observed at low redshift. For example, the 
spectral characteristics of SN 1998ai {z = 0.49) appear to be essentially indistin- 
guishable from those of normal low-redshift SNe la; see Figure 7. In fact, the most 
obviously discrepant spectrum in this figure is the second one from the top, that 
of SN 1994B {z — 0.09); it is intentionally included as a "decoy" that illustrates 
the degree to which even the spectra of nearby, relatively normal SNe la can vary. 
Nevertheless, it is important to note that a dispersion in luminosity (perhaps 0.2 
mag) exists even among the other, more normal SNe la shown in Figure 7; thus, 
our spectra of SN 1998ai and other high-redshift SNe la are not yet sufficiently 
good for independent, precise determinations of luminosity from spectral features 
[73]. Many of them, however, are sufficient for ruling out other supcrnovac types 
(Figure 8), or for identifying gross peculiarities such as those shown by SNe 1991T 
and 1991bg; sec Coil et al. [74]. 

We can help verify that the SNe aX z 0.5 being used for cosmology do not 
belong to a subluminous population of SNe la by examining restframe /-band light 
curves. Normal, nearby SNe la show a pronounced second maximum in the I band 
about a month after the first maximum and typically about 0.5 mag fainter (e.g.. 



[75,18]). Subluminous SNe la, in contrast, do not show this second maximum, but 
rather follow a linear decline or show a muted second maximum [15]. As discussed 

by Ricss et al. [76], some evidence for the second maximum is seen from our existing 
J-band (restframe /-band) data on SN 1999Q {z = 0.46); see Figure 9. However, 
better data on more SNe la are needed to confirm the effect. 
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Figure 7 (left): Spectral comparison (in fx) of SN 1998ai {z = 0.49; Keck spectrum) 
with low-redshift {z < 0.1) SNe la at a similar age (~ 5 days before maximum 
brightness), from [48]. The spectra of the low-redshift SNe la were resampled 
and convolved with Gaussian noise to match the quality of the spectrum of SN 
1998ai. Overall, the agreement in the spectra is excellent, tentatively suggesting 
that distant SNe la are physically similar to nearby SNe la. SN 1994B {z = 0.09) 
differs the most from the others, and was included as a "decoy." 
Figure 8 (right): Heavily smoothed spectra of two high-z SNe (SN 1999ff at z — 
0.455 and SN 1999fv at z = 1.19; quite noisy below ~3500 A) are presented along 
with several low-^ SN la spectra (SNe 1989B, 1992A, and 1981B), a SN lb spectrum 
(SN 1993 J), and a SN Ic spectrum (SN 19941); see [1] for a discussion of spectra 
of various types of SNe. The date of the spectra relative to S-band maximum is 
shown in parentheses after each object's name. Specific features seen in SN 1999ff 
and labeled with a letter are discussed by Coil et al. [74]. This comparison shows 
that the two high- 2; SNe are most likely SNe la. 

Another way of using light curves to test for possible evolution of SNe la is to 

see whether the rise time (from explosion to maximum brightness) is the same for 
high- 2; and low- 2; SNe la; a difference might indicate that the peak luminosities are 



also different [66] . We recently measured the risetime of nearby SNe la, using data 
from KAIT, the Beijing Astronomical Observatory (BAO) SN search, and a few- 
amateur astronomers [77]. Though the exact value of the risetime is a function 
of peak luminosity, for typical low-z SNe la wc find 20.0 ± 0.2 days. We pointed 
out [78] that this differs by 5.8a from the preliminary risetime of 17.5 ± 0.4 days 
reported in conferences by the SCP [79-81]. However, a more thorough analysis of 
the SCP data [82] shows that the high-^ uncertainty of ±0.4 days that the SCP 
originally reported was much too small because it did not account for systematic 
effects. The revised discrepancy with the \ow-z risetime is about 2a or less. Thus, 
the apparent difference in risetimes might be insignificant. Even if the difference 
is real, however, its relevance to the peak luminosity is unclear; the light curves 
may differ only in the first few days after the explosion, and this could be caused 
by small variations in conditions near the outer part of the exploding white dwarf 
that are inconsequential at the peak. 

Extinction 

Our SN la distances have the important advantage of including corrections for 
interstellar extinction occurring in the host galaxy and the Milky Way. Extinction 
corrections based on the relation between SN la colors and luminosity improve dis- 
tance precision for a sample of nearby SNe la that includes objects with substantial 
extinction [28] ; the scatter in the Hubble diagram is much reduced. Moreover, the 
consistency of the measured Hubble flow from SNe la with late-type and early-type 
hosts (sec above) shows that the extinction corrections applied to dusty SNe la at 
low redshift do not alter the expansion rate from its value measured from SNe la 
in low dust environments. 

In practice, our high-redshift SNe la appear to suffer negligible extinction; their 
B — V colors at maximum brightness are normal, suggesting little color excess due 
to reddening. Ricss, Press, & Kirshner [83] found indications that the Galactic 
ratios between selective absorption and color excess are similar for host galaxies in 
the nearby {z < 0.1) Hubble flow. Yet, what if these ratios changed with lookback 
time (e.g., [84])? Could an evolution in dust grain size descending from ancestral 
interstellar "pebbles" at higher rcdshifts cause us to underestimate the extinction? 
Large dust grains would not imprint the reddening signature of typical interstellar 
extinction upon which our corrections rely. 

However, viewing our SNe through such gray interstellar grains would also induce 
a dispersion in the derived distances. Using the results of Hatano, Branch, & 
Deaton [85], Riess et al. [48] estimate that the expected dispersion would be 0.40 
mag if the mean gray extinction were 0.25 mag (the value required to explain the 
measured MLCS distances without a cosmological constant). This is significantly 
larger than the 0.21 mag dispersion observed in the high-redshift MLCS distances. 
Furthermore, most of the observed scatter is already consistent with the estimated 
statistical errors, leaving little to be caused by gray extinction. Nevertheless, if we 



assumed that all of the observed scatter were due to gray extinction, the mean shift 
in the SN la distances would only be 0.05 mag. With the observations presented 
here, we cannot rule out this modest amount of gray interstellar extinction. 

Gray intergalactic extinction could dim the SNe without either telltale reddening 
or dispersion, if all lines of sight to a given redshift had a similar column density of 
absorbing material. The component of the intergalactic medium with such uniform 
coverage corresponds to the gas clouds producing Lyman-a forest absorption at 
low redshifts. These clouds have individual H I column densities less than about 
10^^ cm~^ [86]. However, they display low metallicities, typically less than 10% of 
solar. Gray extinction would require larger dust grains which would need a larger 
mass in heavy elements than typical interstellar grain size distributions to achieve 
a given extinction. It is possible that large dust grains are blown out of galaxies 
by radiation pressure, and are therefore not associated with Lyman-a clouds [87]. 



0.0 i 

^ 0.5 ; 
B 1.0- 



■S 1.5 



SN 1999Q z=0.46 



2.0 r 

l Very Subluminous SN la 
2 5 L (SN 1991bg) 



Typical SN la 
(SN 1995D) 



Subluminous SN la 
(SN 1992bo) 



cq 



1.5 

1.0 
S 0.5 
0.0 
-0.5 

0.50 

„ 0.25 

as 

^^ 0.00 

ca 

^ -0.25 

-0.50 



SN 1999Q 








ml Iiii 


il Normal Dust L 

/// (Av=0.3,nA=0) 

' "Grey" Dust 

(Av=0.3,£2^=0) 

— No Dust 

(Av=0,£J^ > 0) 


























i 




10 


20 30 
Age (days) 


40 


50 



20 40 60 80 

Time (relative to B maximum, days) 



Figure 9 (left): Restframe /-band (observed J-band) light curve of SN 1999Q 
{z — 0.46, 5 solid points; [76]), superposed on the /-band hght curves of several 
nearby SNe la. Subluminous SNe la exhibit a less prominent second maximum 
than do normal SNe la. 

Figure 10 (right): Color excess, Eb-i, for SN 1999Q and different dust models [76]. 
The data are most consistent with no dust and I^a > 0. 



But even the dust postulated by Aguirre [84,87,88] is not completely gray, having 
a size of about 0.1 /im. We can test for such nearly gray dust by observing high- 
redshift SNe la over a wide wavelength range to measure the color excess it would 
introduce. If Ay = 0.25 mag, then E{U - I) and E{B - I) should be 0.12-0.16 
mag [84,87]. If, on the other hand, the 0.25 mag faintness is due to A, then no such 
reddening should be seen. This effect is measurable using proven techniques; so far, 
with just one SN la (SN 1999Q; Figure 10), our results favor the no-dust hypothesis 



to better than 2cr [76], but more work along these hnes is certainly warranted. 

Suppose, though, that for some reason the dust is very gray, or our color mea- 
surements are not sufficiently precise to rule out Aguirre's (or other) dust. If the 
cumulative amount of gray dust along the line of sight grows linearly with increasing 
redshift, we expect that the deviation of the SN la peak apparent magnitude from 
the 1ow-Qm, zero-A model (Figure 3) will continue growing, to first order (Figure 
11). If, on the other hand, the observed faintness of high-^ SNe la is a consequence 
of positive A, the deviation should actually begin to decrease at 2; 0.8 (Figure 
11). In essence, we are looking so far back in time that the A effect becomes small 
compared with Q^m] the Universe is decelerating at that epoch. Thus, we are em- 
barking on a campaign to find and monitor z = 0.8-1.2 SNe la. Given the expected 
uncertainties (Figure 11), a sample of 10-20 SNe la should give a good statistical 
result. 

Note that this test also applies to other systematic effects that grow monoton- 
ically with redshift, as may be expected of possible evolution of the white dwarf 
progenitors (e.g.. [66,67]), or gravitational Icnsing [89]. Indeed, this is our most 
decisive test to distinguish between A and systematic effects. Unless evolution of 
dust, or of the progenitors, or of the lenses is fixed in such a way as to mimic the 
effects of A (e.g., [71]), our claim of > will become much more convincing 
if the deviation of apparent magnitude decreases in the expected manner. Such a 
turnover (Figure 11) can be considered the "smoking gun" of A. 
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Figure 11: The HZT SN la data from Figure 3 (open circles) are plotted relative 
to an empty universe (horizontal line). The two faint curves are the best-fitting 
A model, and the Q^m = 1 (^A = 0) model. The darker curve shows a systematic 
bias that increases linearly with z and is consistent with our z = 0.5 data. The 
expected observational uncertainties of hypothetical SNe la at redshifts of 0.85 and 
1.2 are shown (filled circles). 



CONCLUSIONS 



The luminosity distances of the high-redshift Type la supernovae studied by the 
High-z Supernova Search Team exceed the prediction of a low mass-density (Qm 
~ 0.2) universe by about 0.25 mag. A cosmological explanation is provided by a 
positive cosmological constant at roughly the 3o" confidence level, with the prior 
belief that Qm > 0. We also find that the expansion of the Universe is currently 
accelerating (go < 0, where qo = fln/'^ — ^a)- The independent results of the 
Supernova Cosmology Project are fully consistent with these conclusions. More- 
over, recent precise measurements of the cosmic microwave background radiation 
strongly suggest that the Universe is fiat (l^totai = fin+^A = 1); hence, if f^M ~ 0.3 
(as suggested by many studies, such as of clusters of galaxies), then about 70% of 
the energy density of the Universe must consist of vacuum energy whose precise 
nature and evolution are unknown (but definitely not radiation, normal matter, 
or dark matter). Using the best current values of the Hubble constant, Qm, and 
we find that the dynamical age of the Universe is 14.2 ±1.7 Gyr, including 
systematic uncertainties in the Cepheid distance scale used for the host galaxies of 
three nearby SNe la. This value is comparable to the derived ages of globular star 
clusters. 

Though the consistent results from the microwave background experiments are 
reassuring, we are in the process of testing as exhaustively as possible all systematic 
biases that could be affecting the SN la results. For example, qualitative compar- 
isons of spectra of \ow-z and high-z SNe la do not reveal obvious differences, and 
quantitative tests are in progress. Moreover, the restframe /-band light curves of 
low- 2; SNe la and a single measured high- 2; SN la look similar, as do their broad- 
band colors. The risetimes of low- 2; and high-z SNe la may differ a httle, but the 
statistical significance of this result is not high, and in any case the early part of 
the light curve may have little bearing on the peak luminosity. Further tests are 
in progress. Compelling evidence for acceleration may come in the next few years 
from a comparison of the peak apparent brightness of z ^ 0.8 SNe la with the 
predictions of various models; the signature of nonzero A is quite distinct from that 
of dust, SN evolution, or other effects that grown with redshift. 
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